Diomega production in relativistic heavy ion collisions 
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Using a multiphase transport model, we study the production of a new strange dibaryon (Qfl)o+ 
in dense hadronic matter formed in relativistic heavy ion collisions. The (multi-)strange baryons (E 
and fl) are produced by strangeness-exchange reactions between antikaons and hyperons in the pure 
hadronic phase. The rescattering between the f2s at midrapidity leads to a production probability 
of ~ 3 x 10~ 7 (fi£7) + per event at the RHIC energy of </s — 130A GeV. The production probability 
would be enhanced by two orders of magnitude if (f2fi)o+ and Q, reach chemical equilibrium during 
heavy ion collisions. We further find that the yield of (fifl)o+ increases continuously from SPS to 
the highest RHIC energy. 



PACS numbers: 25.75.-q, 24.10.Lx, 14.20.Pt 

Since first proposed by Jaffe W that the H dibaryon 
with quark content (uuddss) could be a possible light- 
est strangelet — droplet of bound strange quark matter, 
properties of dibaryons have been widely investigated 
in various models, such as the MIT bag model 0], the 
Skyrme model j|||, and the constituent quark model 
p-pf. The consensus regarding the mass of H dibaryon 
is about 2.232 GeV j|. However, extensive experimen- 
tal searches pQ] have not identified any strangelets with 
small values of charge fraction, fz = \Z\/A < 1, and 
strangeness fraction fs — \S\/A < 2. 

Recently, the structure and properties of dibaryons 
with large strangeness are investigated in the chiral SU(3) 
model Jll| [l3[] that have been quite successful in repro- 
ducing several nuclear properties JTJ . In this model the 
H dibaryon is found to be only weakly bound |13j. On 
the other hand, analysis of some six-quark cluster states 
with high strangeness fraction reveals that the diomega 
(f2f2)o+, in particular, is rather deeply bound. Although 
the color magnetic interaction in the one gluon exchange 
term for this system exhibits repulsive feature, the large 
attraction stemming from the chiral quark coupling and 
from the symmetry property of the system lead to a 
rather large binding. In the Resonating Group Method 
calculation, the binding energy of (f2Q)o+ is found to be 
as large as « 116 MeV, and the root mean square dis- 
tance between the two fi's is 0.84 fm. Besides the large 
(negative) charge fraction fz — 1 and strangeness frac- 
tion fg = 3, the new dibaryon (O0)o+ has quite a long 
mean lifetime of ~ 1CV 10 sec as it can undergo only weak 
decay. 

Because of its large strangeness, (f2f2)o+ is not likely 
to be produced in proton-proton collisions. On the other 
hand, strangeness production is enhanced in heavy ion 
collisions and has been suggested as one of the possible 
signals of quark-gluon plasma (QGP) due to large gluon 
density and low energy threshold for ss formation [ ]15| , |l6| |. 
This may thus lead to the formation of exotic deeply 
bound objects composed of quarks or baryons with large 
strangeness. Therefore, in relativistic heavy ion colli- 
sions, especially at RHIC energies, the dibaryon (fifi) + 
could be a new interesting candidate. 



Indeed, recent measurements by the WA97 Collabora- 
tion Q and the NA49 Collaboration Q demonstrated 
substantial enhancement of the (anti-)hyperon yields (A, 
S, and $7) in 158A GeV Pb-Pb central collisions relative 
to p-Pb collisions. The enhancement pattern increases 
with the strangeness content of the (anti-)hyperon. Such 
a large enhancement for multistrange baryons at midra- 
pidity was interpreted as a signal for quark-gluon plasma 
formation. 

On the other hand, even without a phase transi- 
tion (anti)strangeness can be abundantly produced by 
hadron rescatterings alone. In fact, within the micro- 
scopic transport model UrQMD, the WA97 data for mul- 
tistrange baryon enhancement can be explained by re- 
ducing the constituent quark mass in the fragmentation 
of the initial strings in dense matter or by increasing 
the string tension [|l9|. Based on the rate equation ap- 
proach, the multimesonic reactions Y + N <-> wk + nyK 
was demonstrated pOJ to enhance antistrange hyperon 
Y production in a hadronic scenario. Using a multiphase 
transport (AMPT) model, we found that strangeness- 
exchange reactions between antikaons and hyperons in a 
pure hadronic stage also leads to a significant production 
of multistrange baryons at the SPS and RHIC energies 

HI- 

The considerable production of (multi-)strange 
baryons in the dense hadronic stage, in absence of QGP 
formation, could then also result in the formation of 
exotic objects with enhanced strangeness. Properties 
of metastable multistrange baryonic objects consisting 
of nucleons, A's and 5 are studied in the relativistic 
mean- field theory [J22[ . These objects were found to 
have properties quite similar to those of strangelets — 
their quark counterpart. Estimates on the production of 
such strange dibaryons were presented [£3| at the RHIC 
energies by employing wave-function coalescence in the 
RQMD model. 

The pronounced fi production obtained in the AMPT 
model from strangeness-exchange reactions suggests that 
multiple collisions between these omegas in the dense 
hadronic matter may lead to an appreciable production 
of the new strange dibaryon (O0)o+- In this letter, we 



investigate (f2f2)n+ production in relativistic heavy ion 
collisions in the AMPT model. 

The AMPT model || is a hybrid model that uses 
as input both the minijet partons from the hard pro- 
cesses and the strings from the soft processes in the 
HIJING model p5| . The dynamical evolution of par- 
tons are modeled by the ZPC [£6| parton cascade model, 
while the transition from the partonic matter to the 
hadronic matter is based on the Lund string fragmen- 
tation model |P7|| . The final-state hadronic scatterings 
are then modeled by the ART model (28). The AMPT 
model has been very successful in describing the mea- 
sured transverse momentum and rapidity distributions 
of charge particles |29|]3C| l as well as the particle to an- 
tiparticle ratios plfl . The multistrange baryon produc- 
tion is via the strangeness-exchange reactions KK <-> S7r, 
KYi «-> S7r and KB, «-> Cln. Since there is no experi- 
mental information on these cross sections, they are as- 
sumed to be the same as that for KN — ► £71", which is 
the isospin averaged cross section for converting a nu- 
cleon to a sigma |2l| . This assumption is consistent 
with results obtained using SU(3) invariant hadronic La- 
grangians [p2| . For diomega production from collisions 
between omegas, we consider both the the electromag- 
netic process Cl + Cl — ► (fifi)o+ + 7 and the strong inter- 
action process Cl + Cl — > (OQVu. + r\. Their cross sections 
have been evaluated in Ref. |33] | using an effective Hamil- 
tonian. For the process Cl + Cl — > (CICI)q + + 7, the max- 
imum cross section is found to be rs 1.6/ib at pa ~ 0.5 
GeV, while that for the process Cl + Cl — > (17$7)o+ + rj, 
which requires a threshold pn > 0.88 GeV due to finite -q 
mass, reaches a maximum value of ~ 5.5/ib at pa ss 1.1 
GeV. Because of the mismatch between the large initial 
and small final relative momenta between the two omegas 
in the reaction Cl + Cl — > (J7J7)o+ + 7], which involves only 
strong interactions, its cross section is only slightly larger 
than that for the reaction Cl + Cl — > (5757)o+ + 7, which 
involves a much weaker electromagnetic interaction. 

To include the diomega production reactions in the 
AMPT transport model, we parametrize their cross sec- 
tions as functions of center of mass energy y/s (in GeV), 
i.e., 



be unaffected by the production of these particles. For 
(S70)o+ production, since its abundance is extremely 
small, we assume that the production/annihilation rate 
of 77(7) is unaffected. Taking a total scattering cross 
section between two Cls to be cnn = 20 mb, the two 
Cl s will c ollide if their distance of closest approach is 
\J oq£iJti. If Pq and Pq are the production probabili- 
ties of the colliding Cls, the probability of the resulting 
diomega is then P(on) 0+ = p n • Pq • ^n+n-^(Qn) 0+ , where 
rn + n_ > (fi£2) 0+ = (o'n + n_ > (nn)+7 + cr o+fi^(or2)+i))/ f7 nn- 
The probabilities of each of the interacting omegas are 
reduced by Pmnw ■ The annihilation cross section of a 
diomega via the reaction (fiQ)o+ + r\ — *• CI + CI is ob- 
tained using detailed balance of Eq. (g). If a diomega 
collides with a eta, the diomega is totally annihilated 
and two omegas are produced with isotropic momentum 
distribution and each having a probability corresponding 
to the diomega production probability of Pmm , ■ Since 
diomega production is treated perturbatively, the final 
results do not depend on the value of anu used in the 
AMPT model. 
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FIG. 1. Time evolution of midrapidity Q~ and (fin)o+ for 
Au+Au collisions at the RHIC energy of yfs = 130A GeV at 
an impact parameter of 6 < 3 fm in the AMPT model. The 
production and absorption rates for the dibaryon (fifi)o+ as 
functions of time is also shown. 



where ^/s — 2toq = 3.345 GeV, and 



cr n+n^(nn)+i) 



= 1 



^"'('-fd)"* 



(2) 



where ^/s = m n + 2mn — B(oo) = 3.777 GeV. 

Due to small number of multistrange baryons in heavy 
ion collisions, the productions of S and CI as well as of 
the diomega are treated perturbatively in the AMPT 
model H], i.e., the collision dynamics is assumed to 



In Fig. y, we show the time evolution of the abun- 
dances of midrapidity Cl~ and (f2fi)o+ in central Au+Au 
collisions at the RHIC energy of y/s = 130 A GeV. As is 
evident from the figure, most of the omegas and (r2£7)o+ 
are produced within 10 fm/c after hadronization when 
the energy density is high. Beyond this time the diomega 
yield saturates to a value of 2.8 x 10~ 7 per event. 

From the study of time evolution of t he p roduction 
and absorption rates of particles, we found [glj that most 
strange hadrons, including the fis, reach chemical equi- 
librium in collisions at the RHIC energy of y/s = 13QA 



GeV. In Fig. [I] we depict the time evolution of the pro- 
duction (solid lines) and absorption rates (dashed lines) 
of (fiO)o+- It is observed that the absorption rate of a 
diomega by colliding with a eta meson is always smaller 
compared to its production rate at all times. This stems 
from additional cross section in the production channel 
due to the electromagnetic process which has no thresh- 
old. We thus infer that the small number of diomega 
(f217)o+ that may be produced in the relativistic colli- 
sions at RHIC do not reach chemical equilibrium. 
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FIG. 2. Rapidity distributions of E - , QT , and (fifi)o+ for 
Au+Au collisions at the RHIC energy of y/s = 130,4 GeV at 
an impact parameter of b < 3 fm in the AMPT model. 
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FIG. 3. Energy dependence of S~, Q~ , and (f2£7)o+ at 
midrapidity \y\ < 0.5 for heavy ion collisions at an impact 
parameter of b < 3 fm in the AMPT model. 



The rapidity distributions of the produced hadrons 
peak at midrapidity ]24j because of increased multiple 
collisions and rescattering at central rapidities compared 
to that at large rapidities. In Fig. |] we show the ra- 
pidity distribution of the multistrange particles S and Cl, 
along with that of (fiO)o+. The rapidity distribution of 
particles with increasing strangeness content gradually 
become narrower as these particles are successively pro- 
duced from their parent that collide more frequently at 
small rapidities where the baryon density is high. 

We have also studied the energy dependence of 
diomega production in relativistic heavy ion collisions. 
In Fig. pi we show the yield of (Slf2)o+ at midrapidity in 
Au+Au collisions at various energies. Also shown in the 
figure are the yields of Cl~ and S~ as functions of ener- 
gies. It is seen that while the 51 ~~ number increases by 
about a factor of three from the SPS energy (y/s = 17 A 
GeV) to the highest RHIC energy (y/5 = 200 A GeV), the 
(f2!7)o+ number increases by about a factor of 4, indicat- 
ing that the dibaryon (fH7)o+ number reveals a much 
faster rate of increase with y/s as compared to Cl~ . 

The AMPT model gives a lower bound on the diomega 
production probability since there might be other pro- 
duction channels, such as fl - +S~ — > (QCI)q + +K°, which 
have been neglected in the present calculation because of 
their unknown cross sections. If the diomega production 
cross section is larger, then diomegas may reach chemical 
equilibrium with omegas in heavy ion collisions. For most 
other particles including multistrange baryons, equilib- 
rium thermal models have been successfully employed 
p4,B5[ to explain experimental data for their yields and 
ratios over a wide range of energies from AGS to RHIC. 
Adopting the statistical model that is based on the grand 
canonical ensemble with complete thermal, chemical, and 
strangeness equilibrium, we have found that the results 
from the AMPT model for the K + /tt+,P/ P , and K~/K + 
ratios, which are about 0.18, 0.65, and 0.89, respectively, 
at midrapidity in Au+Au collisions at energy ^fs = 1ZQA 
GeV, can be approximately described with a temperature 
T ~ 170 MeV, baryon chemical potential /is ~ 37 MeV, 
and strange chemical potential us — 10 MeV. If we as- 
sume that diomegas are also in chemical equilibrium with 
omegas, then the ratio (ClCl) 0+ /Cl~ is 7.4 x 10~ 5 . With 
the omega number of about 0.41 at midrapidity, this leads 
to a diomega production probability of ~ 3.0 x 10~ 5 per 
event, which is two orders of magnitude higher than that 
obtained in our transport model. 

Diomega production in heavy ion collisions can also 
be studied using the coalescence model [ p3| , p6||3^ | based 
on the omega phase space distribution at freeze out as 
obtained in the AMPT model. Assuming a harmonic os- 
cillator wave function J37J, the Wigner density for the 
diomega is pu(r,q) = 8 exp(— r 2 /d 2 — q 2 d 2 ), where 
r = ri — r-2 and q = (pi — P2)/2 are given in the cm. 
system of Cl-Cl, and the value of d = 0.69 fm corresponds 
to a rms radius of 0.84 fm for the diomega. For cen- 
tral Au+Au collisions at ^fs — 130A GeV, the yield of 
(f217)o+ at midrapidity is found to be 2.6 x 10~ 5 and is 



comparable to that from the thermal model. The similar 
results obtained from both the thermal and coalescence 
models are not surprising as it was shown in Ref. | p8| , 
that the two models are equivalent when matter is in 
thermal and chemical equilibrium and the binding en- 
ergy of the composite particle is much smaller than the 
temperature. Since the AMPT model predicts a hadronic 
matter at freeze out that is close to thermal and chem- 
ical equilibrium, the diomega yield from the coalescence 
model thus should be similar to that given by the thermal 
model. 

Our estimate for the production probability of (O0)o+ 
are well within the limits of the present detectors used 
at RHIC energies. Therefore, this exotic object can, 
in principle, be detected in present and future experi- 
ments. The fact that (fifi)o+ has a large strangeness 
content with a binding energy of ~ 116 MeV, it is stable 
against strong hadronic decays and possess weak decays: 

(oo) + -> 71-- + s° + n- and (oo) 0+ -> 7T° + s- + rr. 

Because three-body decay are involved, the final state 
phase-space would be suppressed. In the sudden approx- 
imation, the mean lifetime of (O0)o+ was found |fL2f to 
be about four times longer than the free O lifetime of 
0.822 x 10 -10 sec. Apart from these conventional decay 
modes, the nonmesonic decay (fifi)o+ — > S~ + 0~ is also 
possible; the estimated U^ lifetime of (fiO)o+ for this 
process is twice the free O lifetime. Thus, instead of di- 
rect observation, the (O0)o+ may also be detected in the 
S~0~ invariant mass distribution. The observation of 
which could provide useful information of the unknown 
— interaction strength. Our study thus opens up 
the intriguing possibility of detecting the new dibaryon 
(O0)o+ in heavy ion collisions at the RHIC energies. 
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